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We investigate the role of CP phases within the renormalizable SUSY 5*0(10) GUT with one 10h, 
one 126// one 126.tr and one 210.jj Higgs representations and type II seesaw dominating the neutrino 
mass matrix. This framework is non trivially predictive in the fermionic sector and connects in a 
natural way the GUT unification of b and r Yukawa couplings with the bi-large mixing scenario for 
neutrinos. On the other hand, existing numerical analysis claim that consistency with quark and 
charged lepton data prevents the minimal setup from reproducing the observed CP violation via 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix. We re-examine the issue and find by inspection 
of the fermion mass sum rules and a detailed numerical scan that, even though the CKM phase 
preferentially takes values in the second quadrant, the agreement of the minimal model with the 
data is actually obtained in a non negligible fraction of the parameter space. We then consider a 
recently proposed renormalizable extension of the minimal model, obtained by adding one chiral 120- 
dimensional Higgs supermultiplet. We show that within such a setup the CKM phase falls naturally 
in the observed range. We emphazise the robust predictivity of both models here considered for 
neutrino parameters that are in the reach of ongoing and future experiments. 

PACS numbers: 12.10.-g, 12.60.Jv, 14.60.Pq, 12.15.Ff 



I. INTRODUCTION 

Testing neutrino properties is one of the greatest chal- 
lenges of the contemporary high energy physics. Thanks 
to the complementary informations provided by solar, 
atmospheric and reactor based neutrino experiments, we 
have obtained in the last few years a consistent descrip- 
tion of the neutrino phenomenology in terms of the phe- 
nomenon of oscillations. The fact that the observed neu- 
trino mass and mixing pattern differs drastically from 
those typical of the quark sector is remarkable and at 
the same time intriguing. On the one hand, the typical 
scale of the neutrino masses (more precisely of neutrino 
mass differences) lies below the electronvolt scale. On 
the other hand, large (almost maximal) mixings appear 
in the lepton sector, at variance with the hierarchical 
structure of the quark sector. The precise data that we 
have obtained in the last years on neutrino properties has 
renewed and spurred the interest on the basic question 
about the origin of fermion masses and mixings, that re- 
mains unanswered within the Standard Model (SM) of 
electroweak interactions. 

The smallness of the neutrino mass may be naturally 
related to a large Majorana mass scale for the right- 
handed (RH) components via the seesaw mechanism Q . 
This mechanism is naturally embedded in grand uni- 
fied scenarios. Of particular interest are Grand Unified 
Theories (GUT) based on 50(10), where all the known 
fermions, including the RH neutrino components, are 
contained in three 16-dimensional fundamental spinorial 
representations. In such a framework the RH Majorana 
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mass is identified with the left-right symmetry breaking 
scale. 

In the last few years a supersymmetric (SUSY) GUT 
model based on the 50(10) gauge group has attracted 
renewed interest for its sharp predictivity of neutrino ob- 
servables [2, b| ■ The model has no more parameters than 
the Minimal Supersymmetric Standard Model (MSSM) 
with massive neutrinos and exact R-parity. In the mini- 
mal renormalizable setting the model contains three gen- 
erations of 16p matter supermultiplcts and the following 
Higgs chiral supermultiplets 0, 0|: 10#, 126#, 126h, 
and 210ir. The 10jy and 126jy representations couple to 
the matter bilinear 16Fi^&Fj — (10s + 120a + 126g)y in 
the superpotential leading to the minimal set of Yukawa 
couplings needed for a realistic fermion mass spectrum 
(5, A denote the symmetry property in the generation 
indices The 126h multiplet also contains a left- 

handed (LH) Higgs triplet that induces small Majorana 
neutrino masses via type II seesaw 0. The 126 h rep- 
resentation is needed in order to preserve supersym- 
metry from Z?-term breaking, while the 210# triggers 
the 50(10) gauge symmetry breaking and provides the 
needed mixings among the Higgs supermultiplets. 

An attractive property of the model is that when dom- 
inance of type II seesaw is invoked the maximality of 
the atmospheric neutrino mixing can be linked to the 
b — t Yukawa coupling unification The model fea- 

tures exact R-parity conservation, due to the even con- 
gruency class (B — L = 2) of the 10 and 126 representa- 
tions (the 120 representation shares the same property), 
with relevant implications for cosmology and proton de- 

cay E El M El EMI m 

In a previous paper [19j we studied in detail the im- 
plications of the mass sum rules of the model on the 
neutrino mass parameters. We showed that in the case 
of CP conserving Yukawa couplings (advocating for the 
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sake of the discussion the soft supersymmetry breaking 
sector as the source of the observed CP violation) , when 
2-cr ranges for the quark masses and mixings are consid- 
ered, the model is consistent with all lepton data at the 
level of the present accuracy. On the other hand, when 
the analysis is limited to 1-er ranges the model shows clear 
tensions in reproducing the observed leptonic spectrum. 
In particular, the electron mass is reproduced by an ex- 
treme fine tuning among quark parameters, while the 
solar mixing angle 9±2 is predicted too close to maximal, 
and the deviation of the atmospheric mixing #23 from 
maximal is too large. The shortcomings of the model 
were previously emphasized in the literature pol l2l| , al- 
beit considering a too limited range of uncertainties in 
the input data. 

In the same paper we worked out a simple renormaliz- 
able extension of the SUSY 5*0(10) model by including 
one 120# supermultiplet to complete the allowed Yukawa 
interactions. On the basis of analytic arguments and nu- 
merical analysis we showed that a small 120# contribu- 
tion to the Yukawa potential allows for a substantial im- 
provement on the fit of the lepton mixings even at the 
l-<7 level. At the same time, the set of the model param- 
eters remains overconstrained by the input data, and the 
framework provides non-trivial outcomes in the neutrino 
sector. 

When complex Yukawa couplings are taken into ac- 
count, recent analyses [2^, I23L |24| show that the fit 
of the lepton mixings may improve, but the Cabibbo- 
Kobayashi-Maskawa (CKM) phase is forced to the second 
or third quadrant by the electron mass fit, thus requir- 
ing significant contributions to CP violation from other 
sources. Non-renormalizable operators are invoked in 
ref. [2^, while in ref. [24| a 120 dimensional Higgs exten- 
sion is considered (with an additional parity symmetry) , 
in order to restore the agreement with the data. CP vio- 
lation in a similar SUSY 5*0(10) framework with a U(2) 
family symmetry was discussed in ref. pl| . 

In this paper we reconsider the study of the fermion 
mass and mixing patterns in the minimal renormalizable 
SUSY 5O(10) model with complex Yukawa couplings. 
We comment on the simple analytic argument that jus- 
tifies the CKM phase taking preferentially values in the 
second or third quadrant and show, by a careful numeri- 
cal treatment of the fermion mass rules, that the minimal 
model exhibit no dramatic tension among quark, charged 
lcptons data and CKM CP violation. The fit of neutrino 
data shows some tension with the value of the strange 
quark mass which is required to be large in order to re- 
produce the recently increased lower limit on the solar 
neutrino mixing. Nevertheless the minimal framework 
does remain a viable GUT candidate 1 . In particular, we 



emphasize the sharpness and the robustness of the pre- 
diction of the U e 3 lepton mixing, that is bound to be 
non-vanishing and within the reach of planned neutrino 
experiments. 

We conclude the discussion by considering the renor- 
malizable 120# extension of the minimal setting pro- 
posed in ref. 01 • While the authors of ref. [24[ , in order 
to sensibly reduce the number of parameters, impose an 
additional parity symmetry that makes all Yukawa cou- 
plings hermitian, in ref. 19] it is assumed that the 120# 
induced contributions to fermion masses are a perturba- 
tion of the minimal scenario (from percent to 0.1 percent 
level). It was shown that the additional (antisymmetric) 
Yukawa interaction leads to a dramatic improvement on 
the fit of the fermion masses and mixings, while main- 
taining, as a perturbation, the predictivity of the mini- 
mal framework. Here we address CKM CP violation and 
show that in the extended model the agreement with the 
data is easily obtained, while maintaining a robust non 
vanishing lower bound for the U e 3 entry of the lepton 
mixing. 



II. FERMION MASSES AND MIXING IN THE 
MINIMAL RENORMALIZABLE MODEL 

Henceforth, unless otherwise stated, we will follow and 
refer to the notation of ref. . When type II seesaw is 
considered, the mass matrices for the SM fermions in the 
minimal renormalizable SUSY 5O(10) scenario are given 
byiEl 

M u = Y 10 v™ + Yuevl 26 , 

M d = Y w vf + Y 126 v d 26 , (1) 

Mi = Y w vf - 3Y 126 v d 26 , 

M v = Y 126 4 26 , 

where Y10 and Y126 are complex symmetric 3x3 matri- 
ces, f.^ 126 denote the VEVs of the components of the 
10jy and 126# multiplets that enter the light Higgs dou- 
blets of the MSSM, and v]? 6 is the tiny induced VEV of 
the LH triplet component in 126 h- These relations can 
be translated into the following GUT scale sum rules for 
the charged lepton and (Majorana) neutrino mass matri- 
ces [12: 

kMi = M u + rM d M v oc M t - M d (2) 

where k and r are in general complex 0(1) functions of 
v u d , and the tilded matrices are normalized to their 
maximal eigenvalue. Let us first rescale by a global phase 



1 Due to the large chiral super multiplets present in the model the 
running gauge coupling diverges shortly above the GUT scale. 
This may call for additional unknown physics to enter below the 



Planck scale |2(I or for an effective gravitational scale to enter 
near the GUT scale. Such a discussion is beyond the scope of 
this paper in as much as perturbativity is maintained up to the 
GUT scale. 
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eq. (J2J such that r becomes real (r = — |r|e l ^ r ) and 
define kl = ke~ z ^ r . Since all the mass matrices are 
symmetric they can be diagonalized by means of uni- 
tary transformations M x — U X D X U X . Taking into ac- 
count that UZU* d = V§ KM = PuVckmP*, where P u = 
diag(e3 iQ i,e^" 2 ,e5 M 3) an d p d = diag(ei i/3l ,e5 i ^, 1) 
parametrize the re-phasing of the quark mass matrices 
necessary to bring the CKM matrix in the standard form 
(denoted by Vckm with one CP- violating phase, Sckm) 
and lead to positive eigenvalues, we may rewrite the sum- 
rules in the M^-diagonal basis as follows 



k' Mi — Vckm D u Vq KM 

k'K cc vS km t d u vS km 



\r\D d 
\r\Dd 



(3) 



,,m& 



The primed matrices in eq. (J3J are given by M' x = 
UjMxUt, while ui accounts for the phase of k' and the 
sign of mf,/m r . The factors h in definitions of phases of 
P u _d are chosen to maintain compatibility with the nota- 
tion used in ref. p^ |. 

One can now exploit the information contained in 
eq. by plugging into the charged lepton mass for- 
mula all the parameters that are known (by running 
the data up to the GUT scale), namely the ratios of 
quark masses in D u = diag(|m u /m t |, \m c /m t \, 1), D d = 
diag(|md/mt|, \m s /rrib\,l), the CKM mixing angles and 
the CP-violating phase Sckm and vary them within 
their experimental ranges. The remaining 8 real pa- 
rameters r|, | A;' |, Arg{k'), at, f3j, that appear in the 
first mass sum-rule, in principle arbitrary are varied over 
their allowed domains. In spite of the many parameters, 
the charged lepton mass sum rule is overconstrained 
and some fine tuning is needed to reproduce the lighter 
eigenvalues. Having obtained an allowed solution of the 
charged lepton mass sum rule, the neutrino masses and 
mixings are then sensitive to the sign of nib /n^r and the 
phase of k' (u>). Since a negative interference between 
the 3-3 entries of the Mi and M d matrices is needed 
to obtain a large atmospheric mixing angle, the phase 
ui becomes strongly correlated to quark phases. As a 
consequence this framework is highly constrained in the 
neutrino sector as well and determines characteristic cor- 
relations among the neutrino parameters. 

In the light of the above remarks it is undoubtedly 
intriguing that this scenario is shown to fit the main 
features of the SM data even when Yukawa phases are 
neglected 0, |2(j. As we have already mentioned, in 
such a case, when considering 1-a uncertainties in the 
input parameters, a rather stringent lower bound for the 
solar mixing angle sin 2 2#i2 > 0.85 appears and the at- 
mospheric mixing angle can hardly be maximal. In ad- 
dition one finds a lower bound for the \U e ^\ parameter 
\U e z\ ^ 0.15, just at the value of the present 90% C.L. 
experimental upper bound \U e z\ < 0.15. While the con- 
straints on the solar and atmospheric mixings are sub- 
stantially relaxed when considering 2-er uncertainties in 
the input data, the lower bound on \U e ^\ is robust and 



represents a characteristic feature of the model as dis- 
cussed in [HI23. 

Before proceeding to the discussion of numerical results 
let us review some issues related to the fit of the fermion 
spectrum in the case of complex Yukawa couplings. 



A. Sckm and the electron mass formula 

When complex Yukawa couplings are considered, it is 
claimed that a successful fit of the electron mass forces 
the CKM phase to take values in the second or third 
quadrant [22I 12H. |2^ | , thus requiring an extension of the 
model to recover the measured CP violation. 

The argument that supports the numerical outcome is 
based on the approximate formula for the electron mass 
eigenvalue that can be obtained from Q3J using the hi- 
erarchical properties of the quark mass matrices in the 
right-hand side (RHS): 



\tfme\e** 1 = -\r\e i01 F d \ A + ■ 
-A 2 A 2 c M3 A 6 ' 
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0(\ 7 



(4) 



Here A = (1 — p — in) where p and n are the Wolfen- 
stein CKM parameters, while F d = ^/A 4 , F c = ^/A 4 
are 0(1) factors. Fitting the normalized electron mass 
(with a typical magnitude of the order of C(A 5 )) amounts 
to compensating the dominant A term in the RHS by 
other terms therein, the only possible one at the given 
order of expansion being that proportional to A 2 . In turn 
this amounts to constraining the size of the CKM phase. 
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The CKM phase is encoded in p and rj as e 
AX 3 (p + in). The typical values of p and n for Sckm 
in the physical region are centered around p ~ 0.21 and 
77 ~ 0.34. Since the parameter lAj 2 = (1 — p) 2 + rj 2 is 
maximized for p < (in refs. [T9L l2p| | the CP conserving 
case n = was considered), the fit of the electron mass 
in formula (0J seems to strongly disfavor the CKM phase 
in the first quadrant |22l| . 

On the other hand, one should be careful in claiming 
the relevance of "subleading" terms in a truncated expan- 
sion. A detailed inspection of the 0(X 7 ) terms in eq. (J3J 
shows that A is not a faithful expansion parameter, in 
that some cofactors, not necessarily dependent on A, can 
become accidentally large (a small denominator in the A 6 
term is an example) . Therefore a larger number of "sub- 
leading" terms in eq. may contribute on top of the 
0(A 6 ) term, and the scan over the complex phases must 
be very detailed not to miss such a solutions (we stress 
that in the minimal model the electron mass is never- 
theless the outcome of a fine tuning; what we are here 
discussing is the extent). As an example let us consider 
the 0(\ 7 ) term 



0(X 7 ) 



A 4 A 2 



J{2a 3 -0 2 ) 



\ r \y 



(5) 
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where F s = ^-/A 3 . Since for typical values of r one 
obtains (1 — \r\) 2 « A 2 , the denominator can be small 
enough to lead to an important correction to the C(A 6 ) 
term. Notice that a small m s favors as well the needed de- 
structive interference in eq. as emphasized in ref. |23j . 

In our numerical analysis we pay particular attention 
to the quality of the parameter scan in those regions 
that may lead to departures from the expectations based 
on the size of A in eq. Once an approximate so- 

lution of the charged lepton mass sum-rule is found (at 
the few percent level) a detailed analysis is performed 
in the neighbor parameter space by linearizing the mass 
relations. Such a procedure improves dramatically the 
convergence of the numerical code, revealing solutions 
that would escape the original scan, unless one performs 
it with extremely high granularity and huge demand of 
computing power. Indeed, such an improved numerical 
analysis shows many solutions of the charged lepton mass 
formula emerging in parts of the parameter space where 
the cancelation among the leading terms in eq. @ is not 
as effective. We discuss and show our numerical results 
in the next subsection. 

The sensitivity of the electron mass fit to sublead- 
ing terms in eq. Q is crucial in the case of the ex- 
tended model considered in ref. with an additional 
120-dimensional chiral super- multiplet. As we shall see 
in Section lllll even for typical magnitudes of the 120 H 
Yukawa contributions to fermion masses as low as 0.1% 
of those coming from the 10 and 126 dimensional Higgs 
multiplets, the role of the A 2 term in eq. J2J is easily 
screened by the 120// induced terms, thus lifting the bias 
on the CKM phase. 



B. Numerical results 

In the first stage we performed a detailed fit of the re- 
lations (|3J). We use GUT scale input data for quark and 
leptons as derived in ref. |28| via two-loop renormaliza- 
tion group analysis. Tablc^shows a sample of quark and 
charged lepton parameters and their uncertainties, for a 
given choice of the supersymmetric threshold scale Ms 
and tan/3. A supersymmetric threshold scale Ms — 1 
TeV and tan/3 = 10, 55 will be considered as typical 
values for the MSSM gauge coupling unification under- 
lying the one-step S'O(IO) breaking scenario here consid- 
ered. The GUT relations between fermion masses and 
Yukawa couplings depends on the vacuum expectation 
values (VEV). Henceforth we follow the MS prescription 
adopted in ref. 28], to which we refer the reader for de- 
tails. 

We consider 90% C.L. ranges for all inp ut parameters. 
With respect to the data used in refs. |28j,|29j we account 
for larger uncertainties in the masses of the light u, d 
quarks and update the range of m s . While m u plays a 
subleading role in the mass sum rules, a light favors 
the reproduction of the electron mass eigenvalue, as we 
discussed in the previous subsection. The values of m u 



TABLE I: Sample of GUT-scale values for quark masses and 
mixings and charged lepton masses. The 1-er data are shown 
for a GUT scale Mq = 2 x 10 16 GeV, a supersymmetric 
threshold scale M S = 1 TeV and tan/3 = 10 pi I2I |30|. 
We use the results of ref. |2H | up to the u, d and s quark 
masses whose central values and uncertainties are updated 
according to the discussion in the text. In our analysis we 
scan over 90% C.L. ranges (1.64 a) of the input parameters. 





0.3585 MeV 


m M 


75.67±°;* MeV 


m T 


1292. 2±^ MeV 


m d 


1.5tg;« MeV 


m s 


29.941^54 MeV 


m b 


1.06t°;S* GeV 


m-u 


0.72±°;^ MeV 


m c 


210.32+2?;a2 MeV 


m t 


82.43±s°;26 GeV 


sin 0i2 


0.2243 ± 0.0016 


sin 023 


0.0351 ± 0.0013 


sin 0i3 


0.0032 ± 0.0005 


ScKM 


60° ± 14° 



and md given in Table I correspond to 4.88 ± 0.57 MeV 
and 9.81 ± 0.65 MeV respectively at 1 GeV 29]. The 
present ranges for the MS running masses at 1 GeV of 
the up and down quarks are 2 to 5.4 MeV and 5.4 to 
10.8 MeV respectively 30]. Accordingly, we will allow 
for values of md at the GUT scale as low as 0.7 MeV. As 
for the strange quark mass an up to date range, which 
includes the lattice evaluations, is m s (2 GeV) = 110 ±20 
MeV, corresponding at the GUT scale to 23 ± 6 MeV. 

The complex phases ^2,3, /9i,2 are sampled in the whole 
range [0, 2tt). The phase u> shows then a correlation to 
the quark phases as a consequence of the tight relation 
between large atmospheric and b — r Yukawa unification, 
which involves a partial cancelation among the terms in 
the RHS of the neutrino mass sum rule. Since the re- 



FIG. 1: The relative density spectrum of allowed solutions 
for the charged lepton masses is shown as a function of the 
CKM phase 8. Although there appear a preference for the 
quadrants with p < 0, there exists a significant number of 
solution with 8c km in the physical region. Quantitatively 
the ratio of the solution density in the preferred area (the 
darkest slice) to that in the 1-a range is about 7 to 1. 
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duced up-quark mass rh u is by far the smallest parameter 
in the mass relations, ct\ does not play any relevant role. 

In Fig. Qwe display the density spectrum of the solu- 
tions of the charged lepton mass formula in eq. Q as a 
function of 5c km in the interval [0, n). 

As one sees the relative density of solutions of the 
charged lepton sum rule is far from being negligible even 
when considering 5 c k m in the 1-a range, contrary to the 
claims in refs. [22,|23. To make this statement quantita- 
tive the ratio of the solution density in the most preferred 
area (the darkest slice in the second quadrant) to that in 
the 1-a 5c km range is about 7 to 1. 

As far as neutrino parameters are concerned, the 
solar mixing angle shows no longer the tight lower 
bound present in the CP-conserving case, namely 
sm 2 20 12 | C p=o > 0.85 (for tan/3 = 10) 19]. On the 
other hand, the experimental improvement on the al- 
lowed values for the solar mixing sharpens a tension with 
the strange quark mass. 



FIG. 3: A density plot of \Am,Q/Am\\ is shown as a func- 
tion of sin 2 2#23 in the minimal renormalizable SUSY 50(10) 
model with complex Yukawa couplings for tan/3 = 10. The 
solid contour encloses the experimentally allowed region at 
the 90% C.L. The dark area corresponds to the solutions that 
are consistent with all neutrino data. 
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FIG. 2: A density plot of sin 2#i2 is shown as a function of 
sin 2 2#23 in the minimal renormalizable SUSY S'O(IO) model 
with complex Yukawa couplings for tan /3 = 10. The solid 
contour encloses the experimentally allowed region at the 90% 
C.L. The dark area corresponds to the solutions that are con- 
sistent with all neutrino data. 

sin 2 20i2 




0.85 



0.95 



FIG. 4: A density plot of \U e 3\ is shown as a function of 
sin 2 2023 in the minimal renormalizable SUSY S'O(IO) model 
with complex Yukawa couplings for tan /3 = 10. The solid 
contour encloses the experimentally allowed region at the 90% 
C.L. The dark area corresponds to the solutions that are con- 
sistent with all neutrino data. 
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0.95 



Fig. shows the area in the solar and atmospheric 
mixing angles that is consistent with the ratio of neutrino 
mass square differences. While the latter (Fig.|3J) bounds 
this area from above, we find that the 90% C.L. lower 
bound on the solar mixing requires a strange quark mass 
above 30 MeV at the GUT scale, that corresponds to 
to s (2 GeV) > 140 MeV (the solutions in the allowed 
region span a GUT scale strange quark mass in the 30—34 
MeV range) . Should the experimental value for the solar 
mixing angle settle above the present central value, it 
would represent a serious shortcoming of the minimal 
S'O(IO) framework. The same conclusion applies to a 
maximal atmospheric mixing. 

The lower bound for the |J7 e 3| parameter is relaxed 
compared to the CP conserving case, although not as 
dramatically as for the solar angle. As Fig. 0] shows, 
for tan/3 = 10 the constraint \U e s\ > 0.15 found in the 
CP conserving setting (see the discussion in ref. [19j) is 



lowered to about \U e 3\ > 0.1. In the case of tan/3 = 55 
we find at the level of accuracy of our numerical analysis 
consistent solutions at the 95% C.L. in the neutrino data. 
On the other hand, the lower bound on \U e 3\ remains 
unaffected. The persistence of such a non-vanishing lower 
bound is a clear signature of the tight correlation between 
lepton and quark Yukawa couplings in this framework, 
that makes |Z7 e3 | ~ 0(X) pf . 

In order to give an explicit numerical example we may 
consider the following (GUT scale) values for the quark 
masses .28] 

m u ~ 0.57 MeV m d ~ 0.73 MeV 
m c - 235.7MeV m s ~ 31.3 MeV 
m t ~ 90.0 GeV m b ~ 1.19 GeV 
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together with the CKM parameters 

sin 12 ~ 0.2253 sinfe ~ 0.0331 
siring - 0.0035 S C km ~ 75° . 



For 





o 144° 




- 142° 


a 3 - 


- 1.2° 




- 0.748 



A' 

U! * 
\k'\ 



216.8° 
224.6° 
-0.2° 
- 0.256 



one obtains the following charged lepton mass matrix 
(normalized to the r mass): 



100 k'M, 



-0.146 + 0.004$ -0.134 + 0.016$ 0.35 + 2.862$ 
-0.134 + 0.016$ -6.985 - 0.197$ 5.321- 11.506$ 
0.35 + 2.862$ 5.321-11.506$ 97.846 + 8.567$ 



TABLE II: Presently allowed values (90% C.L.) for neutrino 
mixing and mass parameters |3l|. 



0.79 < sin 2 2<9 12 < 0.91 
0.93 < sin 2 26» 23 
I sin Sis I <0.15 



2.7 x 10~ 2 < 



< 4.0 x 10" 



lcptons textures. On the other hand, due to a raising 
tension with the strange quark mass a large solar mixing 
and a maximal atmospheric angle can hardly be acco- 
modated. The non- vanishing lower bound for \U e ^\ is a 
robust prediction of the model that falls within the reach 
of the planned long-baseline neutrino experiments (a dis- 
cussion on proton decay and lepton flavor violation pro- 
cesses within this framework is presented in refs. [l7ll23| '). 



The corresponding (GUT-scale) charged lepton masses 
are all within their 90% C.L. ranges: m e = 0.3585 MeV, 
= 75.62 MeV and m T = 1294.0 MeV. The neutrino 
mass matrix is then given by 



100 Ml oc 

' -0.203 + 0.004$ -0.134 + 0.016$ -0.339 + 2.863$ 

-0.134 + 0.016$ -9.404- 0.224$ 5.366- 11.485$ 

v -0.339 + 2.863$ 5.366- 11.485$ 5.855 + 8.952$ 



Keeping into account that the absolute mass scale is set 
by the VEV of the LH triplet in 126/f , one finds neutrino 
mass ratios and mixings (sin 2 2<9 12 = 0.82, sin 2 2# 2 3 = 
0.93, \U e3 \ = 0.11, Am%/Am 2 A = 0.027) that are within 
the present 90% C.L. experimental data. 

The comparison with the data in TablcllTlmust account 
for the running of the parameters from the GUT scale to 
the weak scale. For normal hierarchy with m 2 /m| <C 1 
(which is generally the case in the setup here considered) 
the effects of running of the neutrino mixings and the 
Am 2 ratios down to the weak scale are mild [33|. In 
particular, |?7 e3 | is very stable, with corrections below 1% 
for the whole range of tan j3 considered. The variation of 
the atmospheric angle is small as well, remaining below 
one percent for tan j3 — 10 and at the percent level for 
tan/3 = 55. The largest corrections appear for the solar 
angle at large tan/3: several percents for tan/? = 55. 
Both solar and atmospheric mixing angles grow when 
approaching the weak scale. 

As far as the leptonic CP phases are concerned, the 
Dirac phase Spmns turns out to be generally small 
(< 15°), while the two neutrino Majorana phases y>i,2 
show an approximate 180° correlation. In the example 
reported one finds Spmns = 4°, fi = 10°, ip2 — 191°. 

We conclude that the minimal SUSY 50(10) GUT, 
when complex Yukawa couplings are taken in their gen- 
erality, is not ruled out by present data on the quark and 



III. THE MODEL WITH QUASI-DECOUPLED 
120-DIMENSIONAL HIGGS REPRESENTATION 



In the second part of this paper we consider the case 
of complex Yukawa couplings in a simple extension of 
the minimal renormalizable 5*0(10) proposed in ref. |19| . 
The minimal setting is enlarged to include a 120# chiral 
super multiplct to which the 16p x 16f matter bilinear 
couples. At variance with ref. [24J no a-priori restrictions 
are imposed on the form of the Yukawa couplings. On 
the other hand, in ref. |l9j the 120# contributions to the 
fermion masses are assumed to be two to three orders 
of magnitude smaller than those induced by 10# and 
12Qh- This can be seen as a consequence of a partial 
decoupling of the 120h multiplct and/or a small Yukawa 
coupling. In both cases the setup remains stable under 
quantum corrections and it is (technically) natural. It is 
then shown that, due to the different symmetry property 
of the 120jy Yukawa coupling, the related contributions 
to the mass matrices affect the neutrino mixing angle 
predictions in such a way to improve substantially the 
agreement with the data, even at the l-er level. The 
model does not lose its predictivity and in particular the 
sharp prediction for the U e 3 mixing angle remains. 

In this paper we discuss how the outcomes of the non- 
minimal setting change by switching on CP violation in 
the Yukawa sector. For a detailed discussion of the ex- 
tended setup we refer the reader to ref. • We comment 
here on the features that differ from the real Yukawa case. 

Once the left-handed doublets contained in 10# , 126# 
and 120# acquire a vacuum expectation value (VEV), 
the contributions to the quark and lepton mass matrices 
are generalized to 



M u 
M d 



Y 10 vl Q 



Yl20V d 



(6) 
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M„ 



1211 



12(1 



While Y w and I126 are complex symmetric matrices, 
Yian is a complex antisymmetric matrix. As discussed in 
are taken to be suppressed with respect to 



ref. P3, 



,10 



.,120 

,126 



and v x 2b by two to three orders of magnitude (alter- 
natively the Yukawa coupling should provide the required 
suppression) . In such a case the Yi2o-proportional terms 
in JSJ) can be treated as perturbations of the minimal 
model results. This allows us to maintain the success- 
ful leading order features of the model while treating the 
multi-parameter problem via a perturbative approach. 

The generalized sum rules for the charged lepton and 
neutrino mass matrices then read 



kMi = M u 
M v oc Mi 



tMa 



Yx W ei 



Yi2 {k£i - 

M d 



where 



,120 



Si 



u - re d ) 

Y\2Q£d 



(7) 



m t m b 

As before, one can rotate away the phase of r and diago- 
nalize all the quark mass matrices by means of biunitary 
transformations M x = V^D X V X . The sum rules in 
eq. J7J can be then recast as follows: 

yR^yRj) yO 



k'vfM v V d L 



+ Y{ 20 (k'Ei-e 

« k'vfmvf 



-i<p r 



V 

x 120 



r\D d 
- \r\e d ) 
k'ei 



(8) 



Ife'l 



(D d -Y{ 20 e d ) , (9) 



where k' = ke' 1 ^ , Y[, 



20 



V* ] Y 120 Vf and V^Vf 



Vckm = PuVcKAiPd- Since the antisymmetric compo- 
nents in eq. © are very small, the right-handed quark 

mixing matrix W = V* Vf* can be estimated pertur- 
batively l n the complex case the relevant equation 
reads (x =u,d): 



W 



V CKM 



■2(-\e u \Z^VS KM + \s d \VS KM Z' d 



where Z' x are antihermitean matrices (M T = —M*) 
obeying 

/-ot 



KDu- 
Z'dDd 



DuZ'; 



— L V CKM 



J d 



120 



1 120 V CKM — ^ 



(10) 



and £, 



Re(Z' x )ij - 



Eqs. I|1U|I are then solved by 
Be(Ax)ij 



[Dx)ii \F^x)jj 

-lm(A x )ij 



(11) 



(Dx)ii {Dx)jj 

Having set all relevant notation one can perform a numer- 
ical fit for a given set of the additional parameters Y{ 20 , 
\ s u,d\ and <t> v d - We will follow closely the numerical 
analysis of ref. . 



A. Electron mass formula and screening of the 
CKM CP violating phase 

Using the hierarchical structure of the RHS of eq. (JSJ 
one can again expand the magnitude of the normalized 
electron mass in powers of A: 



l * — Tmm + AT] 



1211 



(12) 



Where the symbol Tmm stands for the minimal model 
contribution (the RHS of eq. J2J). The correction coming 
from the additional terms in eq. (jHJ reads 



ATi 



120 



^Fll(Y{ 20 ) 12 ¥\ 5 + O(\ 6 ). (13) 



where F s = f^/A 3 and F £d = e d /X i are 0(1) form fac- 
tors 

Notice that the C(A 5 ) term of AT120 is in general larger 
than the A 2 term on the RHS of eq. (@J and thus the par- 
tial cancellation of the leading C(A 4 ) in eq. I|12|l is more 
easily achieved. As a consequence, one may expect the 
CKM phase not to be biased towards unphysical values, 
as it happens in the minimal setup. 

By inspection of the 0(A 5 ) term in AX120 and the lead- 
ing 0(A 4 ) term in @, a possible way to make these two 
terms interfere destructively is by taking purely imagi- 
nary entries of the Yukawa matrix Y{ 2n , while assuming 
no spontaneous CP violation. This particular form of the 
coupling is actually obtained in ref. |24j via an additional 
parity symmetry that forces all Yukawa interactions to be 
hermitian. We are now ready to discuss the numerical re- 
sults. 



B. Numerical results 

In analogy with the discussion of the CP conserving 
case we present a sample of the numerical outcomes for 
a given set of the 120_y parameters. According to the 
previous discussion we take 



^120 




\e d \ = HT 3 , lej = 1(T 4 , 



= 



and £; = for simplicity. 

In Fig. [3] we display the relative densities of the 
charged lepton mass solutions as a function of the stan- 
dard CKM phase. As expected, the tension driving the 
CKM phase to the second or third quadrant is almost 
completely screened by the new 120/f induced terms, 
making the physical Sckm a natural outcome of the nu- 
merical scan. Correspondingly, we find that the presently 
allowed ranges for the solar and atmospheric mixing can 



8 



FIG. 5: The relative densities of the charged lepton mass 
solutions are shown as a function of the CKM phase S in the 
extended 120h model, for the setup described in the text and 
90% C.L. input data from Table 



•' &CKM ± Iff 




FIG. 6: A typical density plot of \U e 3 \ is shown as a function 
of sin 2 2023 in the renormalizable SUSY 5*0(10) model with 
complex Yukawa couplings and an additional 120^ Yukawa 
term. The solid contour encloses the experimentally allowed 
region at the 90% C.L. 
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lower bound remains for the \U e ^\ mixing angle, as in the 
CP conserving case 0. As shown in Fig. El obtained 
for 90% C.L. input data and for the 120# setup discussed 
above, we find \U e s\ > 0.05. The bound remains stable 
for tan (3 in the 10 to 55 range. 



IV. CONCLUSIONS 

We have re-examined the role of CP phases within 
the minimal renormalizable SUSY S'O(IO) grand uni- 
fied model, paying particular attention to the fit of the 
charged lepton masses and the neutrino data. We have 
shown, against some prejudice present in previous stud- 
ies, that the observed CP violation in the quark sector 
can be fully accounted for by the standard CKM phase, 
without a severe fine tuning on the charged lepton data. 
While the neutrino Am 2 ratio covers all of the present 
90% C.L. range, both solar and atmospheric angles are 
obtained in their lower ranges. Finally we emphasize the 
sharp prediction for \U e ^\ as a characteristic and robust 
signature of the model. 

As far as the renormalizable extension introduced in 
ref. 0| is concerned, we have found that the fine tun- 
ing in the electron mass is dramatically reduced by the 
presence of small contributions coming from the addi- 
tional Yukawa term. The 120jy induced corrections may 
extend the predictions for the solar and atmospheric neu- 
trino mixings to cover the present experimental ranges. 
In spite of the additional parameters a robust \U e ^ \ lower 
bound remains, that characterizes the mass sum rules of 
the minimal renormalizable S'O(IO) setup. 



be fully covered by tuning the additional contributions. 
The analogues of Figs. and do not add more infor- 
mation and we omit them here. 

It is however worth emphasizing that in spite of the ad- 
ditional freedom in the parameter space a non- vanishing 
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